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Abstract 

We consider the Horava-Witten based model with 5-branes situated near the distant orbifold 
plane and with vanishing instanton numbers on the physical plane. This model has a toric fibered 
Calabi-Yau with del Pezzo base dPj which allows three generations with Standard Model gauge 
group at the GUT scale. Previous analysis showed that the quark hierarchy at the electroweak 
scale could be achieved qualitatively without undue fine tuning due to the effects of the 5-branes 
on the Kahler potential. We extend here this analysis to include the leptons. A new mechanism 
is introduced to obtain neutrino masses by assuming massless right handed neutrinos exist in the 
particle spectrum, which allows a cubic holomorphic term to exist in the Kahler metric, IlH2Vr, 
scaled by the IID Planck mass. After transferring this term to the superpotential, this term gives 
rise to neutrino masses of the correct size at the electroweak scale. With natural choices of the 
Yukawa and Kahler matrix entries, it is possible to fit all mass, CKM and MNS experimental data. 
The model predicts /i — > e -|- 7 decay at a rate that should be detectable for much of the SUSY 
parameter space in the next round of experiments. 
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I. INTRODUCTION 



While the Standard Model (SM) has been successful in fitting all current accelerator 
data, the origin of the quark and lepton mass spectrum remains a puzzle requiring further 
understanding. Thus the explanation of the striking hierarchy of masses (e. g. the up to 
top quark mass ratio is mu/mt ~ 10~^) and the hierarchy of elements in the CKM matrix 
all are beyond the scope of the Standard Model. The matter has been further exacerbated 
by the discovery of neutrino masses, since now in addition there is need for an explanation 
of the MNS matrix as well as the origin of the very tiny neutrino masses. A large number 
of suggestions exist in the literature attempting to explain these properties of quarks and 
leptons. One approach, starting perhaps with the work of Georgi and Jarlskog |l | , suggests 
that the fundamental origin of quark and lepton masses is to be found at high energies, i. e. 
the GUT scale, Mq = 3 x 10^^ GeV, and the complexity we see at low energies arises from 
the running of the renormalization group equations (RGEs) down to the electroweak scale. 
This approach, however, has not appeared to be too promising. Thus a general analysis 
of the u and d Yukawa matrices with five zeros at the GUT scale given in P| is shown 
in Table 1. Here A = 0.2 is the Wolfenstein parameter, and the choice of Table 1, when 
evaluated at the electroweak scale does indeed agree approximately with the quark masses 
and CKM matrix. However, to generate the experimental hierarchy one has to have entries 
at the GUT scale of size ~ 10~^, showing that the problem at the GUT scale is very 
much the same as at the electroweak scale. 



Table 1. Approximate Yukawa textures at Mg for Yu and Yp where A=0.2 



i 



Yrj 



^ v^A^ 0^ 
V2X^ VSA^ A2 



Y, 



D 



^ 2X^0 
2X^ 2A3 
^0 1 



String theory represents at present the only model that has been proposed which in 
principle can calculate the Yukawa matrices from first principles. Unfortunately, math- 
ematical tools to explicitly do this have not yet been developed. In spite of this, the 
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general formulation of the Yukawa problem in string theory opens new windows for 
seeing how the quark and lepton hierarchies might naturally have arisen, approaches 
not available in standard SUGRA GUT theory. In particular, the Horava-Witten het- 

n n 

erotic M-Theory [3|, which offers a natural explanation of why grand unification 



can occur at Mg rather than the Planck scale Mp, has hac 
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27| giving rise 



to three generation models with the SM low energy gauge group SU{?>) x SU{2) x U{1). 
In this model, physical space is one of two 10 dimensional (lOD) orbifold planes separated 
by a finite distance in the 11th dimension, the theory obeying /Z2 symmetry in the 11th 
dimension. Six of the 10 dimensions are compactified to a Calabi-Yau (C-Y) threefold, the 
remaining four being Minkowski space. An array of six dimensional 5-branes perpendicular 
to the 11th dimension can exist between the two orbifold planes. While it is not possible 
to make first principle calculations, one can examine whether the general structure of such 
a theory can replicate the SM at low energy. In this connection, it was seen in j^] that 
the general structure of the quark mass matrices can arise without undue fine tuning if the 
5-branes lie close to the distant orbifold plane, and the instanton number of the physical 
orbifold plane, (3^^^ vanished. It was explicitly shown in 2^ that a three generation model 
with = and SM gauge group indeed can exist for a torus fibered Calabi-Yau (with 
two sections) with del Pezzo base dP-j. The quark and CKM matrix were calculated for a 
model of this type in agreement with data, and it was shown also in approximate analytic 
calculations how the mass hierarchies can arise without undue fine tuning due to the general 
structure of the Kahler potential. 



In this paper we extend the analysis of 



28j in two directions. We first include the charged 



lepton mass matrix and obtain the mass hierarchies experimentally seen. We then consider 
neutrino masses. The conventional way for accounting for the very small mass of neutrinos 
is the seesaw mechanism which gives rise to Majorana neutrino masses. We consider 
here, however, a new way of achieving small neutrino masses based on the structure of the 
Kahler potential. This mechanism is different from the seesaw mechanism, and gives rise to 
Dirac masses for the neutrinos. Neutrino masses and the MNS matrix [30'] are calculated 
consistent with the large mixing anglefLMA) analysis of the solar, atmospheric, reactor and 
long baseline neutrino data (e.g. see 3]| for a global analysis in the context of three-neutrino 
oscillations). Because of the appearance of the MNS matrix, lepton fiavor violation processes. 
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which are absent in the SM, will occur in this model. SM contributions here are too small 
to be observed experimentally. However, SUSY contributions are much larger. We have 
studied the /i — e + 7 decay in this model including all possible contributions. For large 
tan/5, the branching ratio for /i ^ e + 7 is close to the current experimental bound [s^, and 



would be accessible to future experiments j33|, |3^ . A summary of some of the above results 
was given in 3^. 

Our paper is organized as follows. In Sec. 2 we give a brief review of M-Theory, and the 



basic results obtained in 



28| for torus fibered Calabi-Yau manifolds. In Sec. 3 we review and 



update the results of |2S| for the quark masses and extend this analysis to the lepton sector. 
In Sec. 4 we introduce the new mechanism to obtain small neutrino masses and calculate 
the masses and mixing angles for this model. In Sec. 5 we present our calculation for the 
^ e + 7 decay. Finally, conclusions are given in Sec. 6. 



II. HORAVA-WITTEN KAHLER POTENTIAL 



The Horava-Witten M-Theory is concerned with 11 dimensional supergravity on an orb- 
ifold Mio X S'^/Z2, where Z2 is reflection of the 11th coordinate. One can think of this 
space as an 11 dimensional space Mn bounded by two 10 dimensional orbifold planes Mio 
at x^^ = and vrp. In the simplest case, Mio is the product space M4 x X where M4 is 
Minkowski space and X is a (compact) C-Y threefold, the physical world living on one of the 
orbifold planes (e.g. x^^ = 0), the other orbifold plane being a "hidden" sector. In addition, 
there may be six dimensional 5-branes lying along x^^ at bulk points x„ with < a;„ < np, 
parallel to the orbifold planes, each with four dimensions spanning M4, the additional two 
dimensions wrapped around a holomorphic curve in the Calabi-Yau space. The construction 
of a consistent theory involves a remarkable set of interlocking constraints due to anomaly 
cancellation, gauge invariance, and local supersymmetry leading naturally to a theory which 
possesses a number of properties appropriate for phenomenology. Thus there must be Eg 
gauge interactions with chiral multiplets on each Mio orbifold plane (SO (32) being excluded) 
which can easily be broken on the physical plane to the SM group by Wilson lines. The lOD 
gauge coupling constant. A, is uniquely determined in terms of the IID Planck mass, 
leading to the result that the fundamental scale of nature, the IID Planck mass, is 0{Mg), 
and explaining why grand unification occurs at Mq rather than the 4D Planck mass (which 
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is a derived quantity). Finally, a consistent theory exists only as a quantum theory (the 
classical theory being inconsistent), something one would hope might be true for any funda- 
mental theory. Much progress has been made in showing what the low energy structure of 
such a theory might be, and models with three generations of quarks and leptons obeying the 
SM gaugeCToup have been constructed. While the details of the construction of the theory 
given in 3|, y] is rather intricate, it is possible to see how the different elements interact to 
produce a physically interesting model and so we first briefly summarize this construction. 



We then give the relevant 
the latter can be found in 



brmulae needed to examine the low energy structure. Details of 

The field content of IID supergravity is the metric gjj, the gravitino ipij, the three form 
CjjK and its field strength Gjjkl- (In lowest order Gjjkl = djCjKL-) The Bose part of 
the Lagrangian is : 

1^ Jm^^ V 2 48 

where the field strengths obey the field equations Gjjkl = 0, and the Bianchi identity 
dGjjKLM = 0. Here n is the IID gravitational constant. The Horava-Witten theory comes 
about as follows. While in a smooth manifold IID supergravity has no anomalies, on an 
orbifold anomalies arise at the fixed points x = and x = up. To cancel these, it is necessary 
to put Yang Mills multiplets on each Mio orbifold plane, and the cancellation occurs only 
if the gauge group on each manifold is (the phenomenologically desirable) E^. To lowest 
order, the Yang Mills Lagrangian on each Mio reads then: 



4 / d^'x^g tr {-FabF^^ + U^^D, 



AX . (2) 



J^YM - 

where A,B = 1, 2 ... 10, and x is the associated gaugino. However, Eq. is not locally 
super symmetric, and one must proceed in the usual fashion to add additional interactions 
and modifications of the transformation laws to achieve local super symmetry. As usual, 
this involves coupling the gravitino to the Yang Mills supercurrent. However, unlike the 
case where the Yang Mills and supergravity multiplets live in the same space, the gravitino 
here lives in the IID bulk, while the Yang Mills multiplet is constrained to live in lOD. For 
this situation, a locally supersymmetric Yang Mills theory cannot be achieved simply by 



adding interactions on the orbifold plane. It turns out that a supersymmetric theory can be 
achieved only by modifying the Bianchi identities to read 

dGuABCD = 87r2v^^S^+V(")5(xi^ - (3) 

where Xq = 0, x^+i = Trp and Xn, n = 1 . . . N are the positions of the five branes, 



j(o,JV+i) = L (^trFAF--tri?A/?^ . (4) 

167r2 V 2 /,ii=o,.p 

and J("), n = l...N are sources from the 5-branes. With Eq.Q, the total supergravity + 
(modified) Yang Mills Lagrangian can be made locally supersymmetric. However, having 
gained supersymmetry, one has lost Yang Mills gauge invariance (!). For while Eq.(3) implies 
that Gabcd is gauge invariant, the corresponding potential Cuab now is not, i.e. under a 
Yang Mills gauge transformation one has 

5CnAB = (eF^B5(x")) + tr [eFAB^ix^^ - vrp))] . (5) 

which implies the C AG AG term of Eq.iPJ is not gauge invariant. Thus the classical theory 
is not gauge invariant, and a consistent classical theory does not exist. However, in the 
quantum theory, there is in addition the lOD Majorana-Weyl anomaly, and due to unique 
features of the Eg group (!) can cancel the loss of gauge invariance of the "Green-Schwarz" 
C AG AG term provided 

A' = 27r(47r/€)'/l (6) 

Thus only a consistent quantum theory can be built, and this quantum theory determines 
the lOD gauge coupling constant in terms of the IID gravitational constant. 

Eq.® leads immediately to interesting phenomenological consequences. For compact- 
ifying Mn on a Calabi-Yau manifold, one has to lowest order for the 4D gauge coupling 
constant and Newton constant 



where V is the Calabi-Yau volume. Setting V^l^ = I/Mq (so that grand unification occurs at 
the compactification scale as required by the LEP data) and using ac = 1/24, one finds that 
the fundamental IID Planck mass is k~ ^ 2Mg and np-^ = 4.7 x lO^^ GeV. Ahernately 
one may say that the IID Planck mass is the fundamental scale and it sets the GUT scale, 
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while the largeness of the 4D Planck mass is due mostly to accidental Air factors arising in 
the analysis. 



We now summarize the basic formulae of |20[ and 28 1 needed to build a phenomenolog- 
ically acceptable theory. The sources J*^"-* of Eq. (jH)) play an important role in building a 
model. Thus if integrated over a set of independent 4 cycles C^i, they define integer charges: 

= [ J(") (8) 



and Eq.Q then implies S/^l"'* = 0. Here jS^'^^ and jS^^'^^^ are the instanton charges on 
the orbifold planes and (n = 1 . . . A^) are the magnetic charges of the 5-branes. The 
existence of non-zero instanton Yang Mills fields with gauge group G on the orbifold plane 
implies that breaks into G x H where H is the remaining symmetry at the GUT scale of 
the physical theory. We chose here G = SU{5) so that H = SU{5Y. 

Chiral matter arises from the components of the Yang Mills multiplet in the 
Calabi-Yau part of the Mio orbifold [2^. Thus labeling the C-Y indices by holomorphic 
(anti-holomorphic) coordinates a{d) = 1, 2, 3, then one can expand e.g. F^i in terms of a 
basis set of functions in the C-Y space (I is a family index and x a representation index), 
the coefficients in the Minkowski space being the scalar components of the chiral multiplets 
C{Ry^ (where R is the representation): 

F^l = V2^ T,p{R){D^CiR)YP. (9) 

R 

In terms of these quantities, one then defines the metric 

'^^ Jx 

and the Yukawa couplings j20| 

\i.jk{Ri, R2, R3) = [ nA uURi) A U^jiR2) A <(i?3)/Sf ^'""^^ (11) 
Jx 

where Q is the covariantly constant (3,0) form, f projects out the gauge singlet parts, and 
V = vV is the volume of the Calabi-Yau space while v is the coordinate volume: 

V= - [ d^x^; v = [ d'^x (12) 
^ Jx Jx 



^ Recently, an alternate choice, G — SU{4) and H — 5*0(10) has been considered by A. E. Faraggi and 
R. S. Garavuso SEl. 
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In addition one defines the 5", T* and 5-brane moduli by 

Re{S) = V; ReT = V'^^^Ra'; ReZ^ = Zn (13) 

where the modulus R is the orbifold radius divided by p and z„ = x„/ vrp. V can be expressed 
in terms of the a* moduli by Via) = ^dijkd^a^a^ where dijk are the Calabi-Yau intersection 
numbers : 

dijk = / uJiAujj A Uk (14) 

n 

Following the techniques of p, the field equations and Bianchi identities in Eq.(j2I) were 
solved in the presence of 5-branes to leading order 0(fi;^/^) [2^ leading to an effective four 
dimensional Lagrangian at compactification scale Mq- We now state the results that were 
obtained. The gauge kinetic functions on the orbifold planes are given by 

(N 
/3f ) + 5^(1 - Z^f^t^ 
n=l 

= S + e^l^pr'^ + j^Zlpt^^ (15) 

where 

The matter Kahler potential, K = ZjjC^C'^ , on the physical orbifold plane at x^^ = has 
the Kahler metric 



2V 

n=0 



(17) 



where 



and 



Kt = - ln[^d,,fe(r + f%T^ + f^){T'' + T'^)] (18) 



f = - {T + f')Gij - ^{T + f^T' + f^)KTk,T{j (19) 
KT^, = T]j = xf-^ (20) 



The Yukawa matrices are 



YijK = 2V2^XuK ^ l-02\jjK (21) 



for ac = 1/24. The Kahler metric on the distant orbifold plane at x^^ = up is given by 
Eq.lfTTj) with Zn ^ {I — Zn)- 
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III. YUKAWA TEXTURES 



The Yukawa couplings are given in Egs- ljlip and 1)211) as integrals over the C-Y space. 
A priori there is no reason to suggest that a hierarchy such as Table 1 should arise and 
one expects that the non-zero entries to be 0(1). Similarly, one expects a priori that the 
non-zero elements of Gjj in Eq. ()10|) be of 0(1). However, a mild hierarchy can develop in 
the Kahler metric of Eq. (|T7)) if the 5-branes all lie close to the distant orbifold plane, i. e. 
dn = 1 — Zn = 0-1, and provided also /3[^= 0. Then the second term will be small compared 
to the first (e = 0.9), and the model of 



28l | assumed that G/j contributes only to the first two 
generations of the u quark and (which appear together in the SU{5) 10 representation) 
but to all generations of dji, while the second term contributes to all generations but is then 
dominant for the third generation of ul, ur, d^. (That a C-Y manifold exits with (3^^^ = 
with three generations and a SM gauge group is non-trivial and was explicitly shown to be 
possible in 28|.) When the Kahler metric was diagonalized to a unit matrix, it was seen that 
this idea was sufficient to generate a satisfactory explanation of the more extreme Yukawa 
hierarchies at the electroweak scale, and we extend this idea here to the lepton sector. Thus 
the Kahler metric has the general form 



1 



0(1) 0(rf2)\ 



(22) 



0(1) 0(1) 0(rf2) 
\0{d^) 0{d^) 0{d^) J 

where F stands for the different matter fields: q = ul, ur, di, I = (z^l, e^) and e = cr and fx 
is given from Eq.()17p to be e~^'^/'^. We assume that G/j has non-zero elements of 0(1) for 
all generations of dR. (For convenience, we've re-scaled the Zf^ entry in Eq. (22) to 1.) The 
hierarchy then arises when one transforms the Z/j to the unit matrix by a unitary matrix 
U and a diagonal scaling matrix S to obtain the canonical matter fields Cp\ 

(23) 



V JT 



where 



diag5(^) = (A-;/^ A-f ). (24) 

and Ai?i, i = 1,2,3 are the eigenvalues of Zfj/fx- A similar transformation is made on the 
Higgs fields contribution to the Kahler potential 



(25) 



with rescaling of Hi^2' 

= rj^^ K2 (26) 



Before making the transformation of Eq. ()23|). The Yukawa contribution to the superpo- 
tential is [2^ 

Wy = e^^^^YijKC^C-'C^ (27) 

where = ln{S + S) + Kt is the moduh contribution to the Kahler potential. From Eqs. 
(fni) and (dHl), one has 

= -ln{2V) - ln{8R^). (28) 
Written in terms of SM fields Wy then is 

= 4R37^(^^"^^i^2Wi? + Y^'^kLH.dn + Y^^hH^en). (29) 

and after the transformation to the canonical matter fields one has 

Wy = M^A^^^M'^iJ^ + d'L\^'^WjiH[ + e^A(^)e^/7(. (30) 

where A'-"''^''^-' are give by ^ 

>-u = ^;jp^-^(Sl«>C/«y<-'C/<"'S<"))„ (31) 

We use here the notation "~" for transpose. In Eq. ()3()|1 . A*-"''^'^-' play the role of the Yukawa 
matrices at the GUT scale in the phenomenological analyses such as in However, in 
general they are not symmetric matrices and so M-Theory textures are uniquely different 
from what has previously been considered in phenomenological analyses. In brief, it is 
the smallness of the third generation eigenvalues of the Kahler matrices appearing in the 
denominators of Eq. ()31H3^ (from the factor S of Eq. ()24|) ) that give rise to the large third 
generation masses. 



^ We correct an error in (23, the omission of the V 3 factor in Ea. (|13|l (see e.g. |38|), which leads to a 
factor 1/1/^5 in Ea. H31l33|l rather then \/V. 
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Table 2. Kahler matrices z'f'j'^'^'^^ and Yukawa matrices for tan/3=40. 



1 0.3452 
0.3452 0.1311 0.006365 
0.006365 0.00344 



1 -0.547 
-0.547 0.432 0.025 
0.025 0.09 



T 



I 1 0.496 ^ 

0.496 0.564 0.435 
y 0.435 0.729 



I 

I 1 0.624 ^ 

0.624 0.397 0.00574 
y 0.00574 0.004407y 



diagy^"*) = (0.0114, 0.0597, 0.104 exp[0.657ri]); 
diagy^"') = (2.052, 0.2565, 1.8297); 
diagy('=) = (0.307, 3.789, 1.821). 

In j^l we saw for the case of tan /3 = 3 how the above Yukawa matrices gave rise to 
the experimental quark masses and CKM matrix elements at the electroweak scale, and we 
showed there analytically how the hierarchies arose naturally without undue fine tuning. 
We now update this analysis for the case of tan /3 = 40, and extend the discussion to include 
the lepton sector. Table 2 shows a choice of Kahler metric and Yukawa matrices that satisfy 
all the current experimental data. The Z^, and entries for F = u,l,e are 0{d?') 
(for d = 0.1) as required by Eq.()22|). For simplicity we have assumed that the q and u 
quarks have identical Kahler matrices and have the maximum number of zero entries, and 
that the Yukawa matrices are diagonal. One phase is assumed in the Yukawa matrices to 
account for CP violation. To compare with low energy data, we use one loop Yukawa RGEs 
and two loop gauge RGEs to evaluate the Yukawa couplings at the electroweak scale, which 
we take to be m^. Below rrit we assume that the Standard Model holds and include in 
our calculations the QCD corrections (which are quite significant). The QCD correction 
factors used were rjc = 2, 1]^ = 2.5 = rid, r]h = 1.6 and r]s = 2.5. Diagonalization of the 
low energy Yukawa matrices then allows one to generate the low energy quark and lepton 
masses and the CKM matrix elements. The results are shown in Table 3, and are in good 
agreement with experiment. Of course in a fundamental analysis, the precise entries in 
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Table 3. Quarks and leptons masses and CKM matrix elements obtained from the model of Table 
2. Masses are in GeV. Experimental values for lepton and quark masses are from [32] and CKM 
entries from y unless otherwise noted. 



Quantity 


Theoretical Value 


Experimental Value 


mt(pole) 


175.2 


174.3 ±5.1 


mc{mc) 


1.27 


1.0-1.4 


mu{l GeV) 


0.00326 


0.002-0.006 


mb{mb) 


4.21 


4.0-4.5 


m,(l GeV) 


0.086 


0.108-0.209 


md(l GeV) 


0.00627 


0.006-0.012 


171^ 


1.78 


1.777 




0.1054 


0.1056 


nie 


0.000512 


0.000511 


\Vus\ 


0.221 


0.2210 ±0.0023 


\Vcb\ 


0.042 


0.0415±0.0011 


\Vub\ 


4.96 X lO"-'^ 


Smtoit ± 0.45 X IQ-^ 


\Vtd\ 


6 X 10-'^ 


9.2 ± 1.4 ±0.5 X 10-^ 


sin 2/? 


0.803 


0.731 ± 0.056 HO] 



Table 2 arise from integrals over the Calabi-Yau space, an analysis that cannot at this stage 
be performed. However, our discussion has shown that the general structure of the Kahler 
metric and Yukawa couplings arising in our Horava-Witten model can lead to low energy 
quark and lepton spectra consistent with all current experiments without the fine tuning 
used in phenomenological analyses. 

Without knowledge of the value of the factors R^V^^'^ ^yGH^ in the denominators of 
Eq. ()31ll3^ . Kahler textures can only determine the mass ratios. As in 2^, we use the top 
Yukawa at the GUT scale to determine the value of this common factor. If we write V = r^, 
where r is the mean radius of the Calabi-Yau manifold divided by the co-ordinate radius. 
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then for 6*^1,2 = I5 finds that 

Rxr = 6.82. (34) 

In the next section, we will show that R and r can be determined separately if massive 
neutrinos enter our model via the mechanism proposed there. 



IV. NEUTRINO MASSES AND OSCILLATIONS 



In the last section we presented a way to generate the Yukawa textures in the quark and 
lepton sectors whose structures are the same as the SM. The consequence of the masslessness 
of neutrinos in the SM is that the mass eigenstates of leptons are identical to their gauge or 
flavor eigenstates and, unlike the quark sector which has a CKM mixing matrix, the lepton 
sector does not. Therefore, there is no oscillations between neutrinos in the SM. However, 
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m, SNO 



43| and KamLAND 
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the neutrino experiments of Super-Kamiokande 
have shown the existence of neutrino oscillations which indicates that neutrinos are actually 
massive particles. In this section we will show that massive neutrinos can be included in 
our model and their masses and mixings can be fitted into the large mixing angle (LMA) 
solution. (For a recent review of neutrino oscillations see |l45l| . ) 

The simplest way to include massive neutrinos to our model is to associate a right-handed 
neutrino to every left-handed neutrino and insert by hand a term proportional to 

Y^'^HlH^pr (35) 

into superpotential (j^HI)- However, the Yukawa couplings in the neutrino sector have to 
be extremely small and thus this solution is theoretically less interesting unless there is a 
mechanism behind it. The most widely used way to overcome this problem is the seesaw 
mechanism 2^. In seesaw models, besides the usual Dirac mass terms (which are ap- 
proximately the same size as other fermion masses), one introduces additional very large 
Majorana masses which enter in the off-diagonal entries of the neutrino mass matrix. As a 
consequence, some eigenvalues are suppressed to the desired values when the diagonalization 
of neutrino mass matrix takes place. The physical neutrinos in seesaw models are then of 
Majorana type while other leptons and quarks are Dirac fermions. Here we propose a new 
way to generate neutrino masses. In our model, neutrinos are of Dirac type and thus the 
similarity between leptons and quarks is preserved and no neutrinoless double beta decay 
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exists. We will see that our new mechanism provides a reasonable physical explanation to 
the origin of term (jH^j) . 

The Kahler potential in principle can have gravitationally coupled trilinear terms which 
are usually ignored as they generally are of negligible size. However, we assume here that 
our Kahler potential at the GUT scale contains the holomorphic cubic term K^^^ = K^ + KJ 
where 

K, = m^Y^''nLH2iyR (36) 

where 1 / nu is the 11 dimensional Planck mass (i.e. 1 / nu ^ Mq) and Y" is a Yukawa matrix. 
We note that Eq. (|36|) is the only gauge invariant holomorphic cubic lepton term involving Upi 
and that nu is the natural scale for Horava-Witten theory. We assume here that the Yukawa 
contribution to the superpotential is still given by ((201), ^^at no additional neutrino 
masses arise there. One can transfer K,^ from the Kahler potential to the superpotential by 
a Kahler transformation {I/k^ is the 4D Planck mass): 

K -> K-K^^\ 

m _^ (.^iK^W = W + K^K^W + ■■■ (37) 

Now when supersymmetry breaks, the superpotential W will grow a VEV of size: 

{W) ^ X^Ms (38) 

where M5 is of electroweak size. Consequently, after supersymmetry breaking, an additional 
term appears in superpotential (jHUj) : 

f K<".;.//..„. (39) 

In the above we have assumed that no additional neutrino masses arise in the superpo- 
tential. One might imagine that this could come about if there exists a global symmetry or 
a non-gauge discrete symmetry that obeys. (Examples might be if was a member or a 
global SU{2)r doublet or were charged under a global U{1) symmetry; or if appeared in 
W only as an even power one could assign it a discrete quantum number —1.) Then in the 
cubic part of the superpotential, the global or discrete symmetry (plus SU (2) l symmetry) 
would forbid Majorana or Dirac masses formed from vl.r- Since gravity is expected to break 
such global or discrete symmetries in string theories, quartic terms scaled by kh could arise 
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(which is in fact why the contribution of Eq. ()36p would be allowed in the Kahler potential 
in this scenario.) Such terms would have the general form {H2IY and HiH2i^% giving rise 
to Majorana masses. Using the parameters of Table 2 and Eq. (|43|) below, we estimate that 
the first term gives rise to neutrino masses a factor of ~ 100 smaller than Eq. (j36|) . and 
the second term is a factor of 1/ tan/? smaller yet. Thus they would produce only small 
corrections in our model. We have not investigated in detail if one could dispense with 
Eq. ()36|) completely and try to get the correct neutrino masses using only the above super- 
potential terms, but because of the above result, we believe it to be unlikely. We note also, 
that the above superpotential terms could also occur in the standard see-saw model, but 
are neglected there. While global and discrete symmetries can indeed arise from Calabi-Yau 
manifolds in Horava-Witted M-theory, whether or not the necessary symmetries are present 
in physically interesting Calabi-Yau manifolds is not known. Thus we have postulated their 
existence in this paper. ^ 

We can now proceed as in Sec. 2. First diagonalize and rescale the Kahler matrices Zu 
of vr and other fields to the unit matrix. Then make the necessary transformations in the 
superpotential to the canonical normalized fields. The term giving rise to neutrino masses 
can then be written as 

',\^^^u'^H', (40) 



where 

Note that the overall coefficient in (PT|) is different from the one in ()31II33|) because the 
neutrino term originates from the Kahler potential, not the superpotential (jTfj) which has 
the additional coefficient e^^™. It is thus possible to use the experimental neutrino mass 
square differences to determine R. In the example given below, we find that R = 2.13 
produces acceptable neutrino masses (we assume M5 = 1 TeV in our calculation), and from 
Eq.()34j). one finds that r = 3.20. At the weak scale, after the diagonalization of charged 
lepton and neutrino Yukawa matrices, the Maki-Nakagawa-Sakata (MNS) lepton mixing 
matrix arises. We follow the standard parameterization Q (the phase similar to the one 



We thank the Referee for bring this point to our attention. 
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in the CKM matrix is ignored): 



V, 



MNS 



C12C13 S12C13 Si3 

~'5i2C23 — C12S23S13 C12C23 — S12S23S13 S23C13 
■512-523 ~ C12C23S13 — C12S23 — S12C23S13 C23C13 



(42) 



where Cij = cos%, Sij = sin 6ij and i,j = 1,2,3. 



The following is an example at tan j3 = 40. We use the lepton entries of Table 2, and the 
following neutrino Kahler and Yukawa matrices at Mq: 



( 



1 -0.465 
-0.465 0.3105 0.0254 
0.0254 0.027 



(43) 



diagF^^) = (4, 0.4, 4). 



(44) 



The neutrino mass square differences and mixing angles at the weak scale are then calculated 
to be: 



= 5.5 X 10^"^) eV^; 
Amj^ = 2.7 X 10^-^^ eV^; 



(45) 
(46) 



tan^ 9i2 = 0.42; tan^ ^23 = 0.93. (47) 

with \Ue3\ = 0.005. Since our model is a complete model of neutrino masses, we can calculate 
all the masses themselves and not just the mass square differences. For the above example 
we find 

mi = 6.5 X 10"^ eV; m2 = 7.4 x 10~^ eV; = 5.2 x 10"^ eV (48) 

consistent with cosmological constraints on neutrino masses {47 1. 

The analysis of solar and KamLAND data in terms of two neutrino oscillations gives for 
the LMA solution m 



0.20 < tan^ Og < 0.68 ; 5.6 x 10"^ < Am|/eV^ < 8.9 x 10" 



(49) 
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where Am| is the solar neutrino mass square difference and 63 is the corresponding mixing 
angle and the ranges in Eq. ()49|) (and Eq. ()50|) below) are 3cr around the central value. The 
analysis of Super-Kamiokande and K2K data shows for the LMA solution 48|: 



0.85 < sin^ < 1 ; 1-4 x 10"^ < AmyeV^ < 3.8 x 10"^ (50) 

where Am^ and 6a are the relevant mass square difference and mixing angle for the atmo- 
spheric neutrino oscillation. 

Since ,„ ou. ca.e \U,, . 0. sola. a.d at„„.pKe* „eutn„o osc.llat.o. decouple HQ. 

Therefore the two neutrino oscillation analysis can be applied to our case with the effective 
mixing angles given by: 

6s = ^12 , 6a = 623- (51) 

Eq. (pSj) - (P7j) . and (plj) show that our results agree with the current LMA solution quite 
well. 



V. /i ^ e + 7 DECAY 

Lepton flavor violation (LFV) processes in supersy mmetric models have been discussed 
in much detail in the literature (e.g. 0, 0, 0, 0, 13)- '^^^ model the MNS matrix 
by itself can give rise to LFV processes, but this contribution is still very small, e.g. for 
yU — > e + 7, the decay rate is of (9((mj,/mH^ Therefore, In this section, we 

mainly discuss the additional supersymmetric contributions to the LFV process /i — e + 7. 

The operator for yU ^ e + 7 is: 

ic 

£^^e+7 = e (^'"'Qu {aiPh + arPR) f^-A^ + h.c. (52) 

where Pl,r = (1 =F75)/2 and a^'^ = | [7^, ■j'^]. The decay width for /i — >• e + 7 can be written 
as: 

r(/i^e + 7) = ^(|a/r + |a.r) (53) 
Then the branching ratio is given by: 

T{fi ^ e + 7) 



Br(yU ^ e + 7) 



«i|' + l«r|') (54) 
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FIG. 1: Feynman diagrams for the neutralino and chargino contributions to ^ e + 7. 

The supersymmetric contributions include the neutrahno and chargino diagrams shown 
in FIGffl The neutrahno diagram gives: 

4 6 _2 



E E ^ i^fL^fL^.^^.) + ^N-*N- F.ixr,)) (55) 
i=i k=i Ik ^ ^ ^ 

+ '-^K-NilF,(.g)) (66) 

i=l k=l '■ ^ ^ 



where xfj. = m -0 / mf and 



Fi(x 



2a;^ + Sx^ — 6a; + 1 — 6x In a; 
12(x - 1)4 
— 1 — 2x\nx 



2{x-lf 

where (71 and g2 are the U(l) and SU(2) gauge couphng constants, X is the matrix diagonal- 
izing the 4x4 neutrahno mass matrix M^o according to M^oX = X*M'^\ D diagonahzes 
the 6x6 charged slepton mass matrix Mf according to M?D = DM^^^ and F^(^) is the 
diagonahzed Yukawa matrix of charged leptons (We use the notation of [59]). 
Similarly, the chargino diagram gives 

2 6 3 o 



a? 



E E E (cfr C£F3(4) + !!!ltcfrc« F.(.g)) (68) 

i=l k=l m=l i^k ^ A* / 

2 6 3 ^2 / ^ ^ N 

E E E C'ik*C',hFsixg) + ^C^lCcllF.ixg) (59) 

i=l k=l m=l '^k ^ / 
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where = m^±/m|^ and 



F3(x) = 



— + 3a; + 2 + 6a; In x 



12(a; - 1)4 
x'^ — Ax + 3 + 2\nx 



Cfj? = Ut2Pmk{^MNs)'lmXi 



iD) 



(60) 



where Vmns is the MNS mixing matrix, P is the matrix diagonahzing the 6x6 sneutrino 
mass matrix M| according to M?P = PM1^^\ U and V diagonahze the chargino mass 
matrix M^± according to U* M.^±V^ — and y^-^) is the diagonahzed Yukawa matrix 

of neutrinos. 

To evaluate the branching ratio of — > e + 7, we first generate Yukawa textures in the 
way described in Sec. 3 with phenomenological inputs including the fermion masses and 
neutrino oscillations described in Sec. 4. Then wc choose soft breaking parameters at the 
GUT scale and run the RGEs to the weak scale. Finally one can use the formula given 
above to calculate the // — > e + 7 branching ratio. We display our results in the following 
three figures for tan /3 = 10, 30 and 40. 




FIG. 2: Branching ratio of /x — > e + 7 for tan/? = 10 

Although our Yukawa textures are constructed through the Kahler potential, the usual 
mSUGRA structure holds at the GUT scale. In addition, to constrain the parameter space. 
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FIG. 3: Branching ratio of /i — > e + 7 for tan /3 = 30 
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FIG. 4: Branching ratio of // ^ e + 7 for tan /? = 40 



we implement the relic density constraint 



47|: 0.095 < n^oh^ < 0.129. For a given mi/2, an 



6OII . Since the 



allowed narrow region of mo is determined by the relic density constraint [59 
/i ^ e + 7 branching ratio is not sensitive to the value of mo we only show the result for one 
value of mo in the allowed region for any given mi/2- There are other experimental bounds 
on the parameter space, e.g. the 6 — s + 7 decay, the light Higgs mass, the muon g — 2, all of 
which can easily be implemented. Since mo is significantly constrained by the relic density 
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bounds for any given mi/2, other experimental constraints are only needed for constraining 
mi/2- For example, at tan/5 = 40, the 6 s + 7 branching ratio produces a lower bound of 
~ 400 GeV on mi/2 while a muon g — 2 deviation from the SM can produce an upper bound 
on mi/2- For the purpose of showing what yU ^ e + 7 branching ratio can be reached in our 
model, except for the relic density constraint, we ignore the other experimental bounds in 
our plots since they are not significant for this purpose. 

In our plots the y-axis is the logarithmic ratio of our theoretical predictions to the current 
experimental bound js^. Therefore, only the region below zero is experimentally allowed. 
One can see that at large tan/3, especially for tan/3 = 40, the theory predictions are only 
about one order of rnagnitude smaller than the experimental bound and hence accessible to 
future experiments js^ 13 while part of the parameter space for lower tan (3 will also be 
accessible. 



VI. CONCLUSION 



In this paper we have extended a model of the quark mass hierarchy based on the 



Horava-Witten M- Theory 28| to include charged leptons and massive neutrinos. The model 
is based on the assumptions that five branes exists in the bulk lying near the distant orbifold 
plane (i. e. about 90% of the way from the physical plane), and that the instanton charges 
on the physical plane vanish. We had previously seen that this gave rise to a three gener- 
ation model with the Standard Model gauge group at the GUT scale. While one cannot 
calculate Yukawa couplings in M- Theory (they involve integrals over the Calabi-Yau space) 
these constraints were sufficient to qualitatively account for the quark mass hierarchy at the 
electroweak scale without undue fine tuning. The mechanism that achieved this was that the 
five brane contribution to the Kahler potential gave rise to small Kahler matrix eigenvalues, 
and the quark masses were proportional to the reciprocal square root of the eigenvalues 
when the kinetic energy was put into canonical form. We saw that the same mechanism 
also gave rise qualitatively to the hierarchy of charged lepton masses, again without any 
excessive fine tuning. 

Neutrino masses can arise in these models if a right handed neutrino exists in the massless 
particle spectrum. Then one can assume that the Kahler potential has a cubic holomorphic 
contribution of the form of Eq. (p¥)|) . the interaction being scaled by the 11 dimensional 
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Planck mass (the basic parameter of Horava-Witten theory). When transformed to the 
superpotential by a Kahler transformation, this term gives rise to neutrino masses of the 
correct size after supersymmetry breaking. (Thus the mechanism being used here for the 
neutrino masses is similar to the one previously used to generate a fi parameter of electroweak 
size 0]-) it is possible then to chose natural sized values for the Yukawa and Kahler matrix 
entries to generate masses and CKM and MNS mixing angles in agreement with all low 
energy data. The neutrinos in this model are Dirac, and so will exclude neutrinoless double 
beta decay. However, the mixing in the neutrino sector allows for /i — e + 7 decay to occur, 
and with reasonable values of the SUSY parameters, this decay should become observable 
in the next round of fi catalysis experiments [ssl . 3^ over a significant range of parameters. 

Aside from the Kahler and Yukawa matrices, the quark, lepton and neutrino properties 
depend on the Calabi-Yau volume modulus V which we have parameterized by V^^^ = r and 
the radius modulus R. We have found that all the quark, lepton and neutrino masses can 
be fit satisfactorily with r and R of 0(1). Thus for the example in text for tan/5 = 40 we 
found R = 2.13 and r = 3.20. One important feature of this Horava-Witten model that has 
not been addressed here is how to stabilize the position of the 5-brane close to the distant 
orbifold plane. One possibility may involve quantum corrections, e.g. membrane potentials 



between the 5-brane and the orbifold planes 
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